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Thomas  S.  Mautncr 
Naval  Ocean  Systems  Center 
San  Diego.  CA  92152 


ABSTRACT 


Based  on  the  work  of  Greenblatt  (1978).  an  enhanced  optimiza¬ 
tion  technique  for  the  reduction  of  propeller  unsteady  forces  and 
the  determination  of  skew  distributions  has  been  developed.  The 
current  method  provides  an  efficient  propeller  design  tool  capa¬ 
ble  of  determining  a  variety  of  cubic  or  quadratic  skew  distribu¬ 
tions.  subject  to  constraints,  which  minimize  the  unsteady  forces 
produced  by  the  various  harmonic  components  of  the  input 
wake.  The  original  skew  optimization  method  was  extended  to 
include  higher  order  harmonics,  and  the  original  force  calculation 
method  was  replaced  by  an  expended  version  of  the  method 
developed  by  Thompson  (1975).  Calculation  of  forces  and  skew 
distributions  associated  with  a  representative  propeller  show  that 
acceptable  reduction  of  unsteady  forces  can  be  obtained  without 
having  to  place  severe  restrictions  on  the  model  constraints. 
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INTRODUCTION 

One  current  and  important  issue  in  the  design  of  both 
marine  and  aircraft  propellers  is  the  reduction  of  propulsor  gen¬ 
erated  noacTHl  men*  pngUcr  designers  ate  concerned  with 
die  reduction  of  noise  due  to  both  the  transmission  of  unsteady 
propeller  forces  through  shafting  into  the  vehicle  and  propeller 
noise  radiated  into  die  near  and  far  fields.  For  example,  the 
reduction  of  self-noise  in  a  marine  vehicle  will  result  in  improved 
sensor  operation  while  the  reduction  of  aircraft  noise  improves 
cabin  comfort  and  reduces  ground  noise. 

- The  problem  considered  herein  is  the  generation  of 

propeller  noise  due  to  non-uniform  inflow  velocity  fields.  This 
has  been  and  still  is  an  issue  of  research  in  the  area  of  marine 
propellers,  and  it  is  felt  that  a  recently  developed  method  of  cal¬ 
culating  unsteady  forces  -and  xoccilVint  skew  distributions  (Le. 
blade  sweep)  for  marine  propehcnr’bny  also  be  applicable  in  the 
design  and  analysis  of  aircraft  propeilers.xin  particular  pusher 
propellers. 'Although  the  sources  of  the  .flow  non-uniformities 
may  be  somewhat  different  for  marine  And  aircraft  propellers, 
flow  non-uniformities  will  result  in  tbe  generation  of  unsteady 
forces  and  noise  irregardlcss  of  the  propeller  type.  For  example,  a 
propeller  operating  on  an  axisymsnbtric  body  with  appendages  (a 
typical  underwater  case)  experiences  a  non-uniform  inflow  velo¬ 
city  field  due  to  the  upsteam  appendages  (for  example  Nelson 
and  Fogarty,  1977)  while  uf  aircraft  prop-fan  is  subjected  to  flow 
non-uniformities  due  to  the  wing  installation,  the  nacelle  or 
engine  exhaust  and  a  busker  propeller  is  exposed  to  s  non- 
uniform  flow  field  dub  so  the  upsteam  airframe  geometry  (for 
examples  see  Meogbr,  19*4;  Schulten.  1914;  Meager  and 
Rohrbach,  19*6;  and  Takallu  and  Mock  19*7). 
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The  generation  or  unsteady  propeller  tllade  forces  can 
be  described  by  consideration  of  a  propeller  operating  in  the  tur¬ 
bulent  wake  of  a  body  having  upsteam  protuberances  which 
create  wake  non-uniformities  which  in  turn  result  in^patial  and 
temporal  fluctuations  of  blade  angle-of-anack.  Thcsc*angic-of- 
attack  fluctuations  result  in  unsteady  blade  loadings  and  the  gen¬ 
eration  of  propeller  noise,  and  the  noise  sources  are  character¬ 
ized  by  three  types  of  unsteady  force  mechanisms-  a)  turbulence 
injestion;  b)  vortex  shedding;  and  c)  Wadc-rate^Thc  first  two 
mechanisms  typically  generate  continuous  spcctrum\broadband) 
radiated  noise  while  blade-rate  forces  generate  discrete  fre¬ 
quency  noise  at  various  blade-passage  frequencies  and\harmon- 
ics.  This  paper  will  address  the  reduction  of  blade-rat^  noise 
which  provided  the  original  motivation  for  the  application  of 
skew  in  propeller  design.  \ 

An  example  of  a  non-uniform  wake  incident  upon  a 
propeller  is  one  generated  by  an  axisymmetric  body  with  append¬ 
ages.  The  boundary  layer  behind  the  appendage  (i.e.  figure  1)  is 
characterized  by  a  complex  velocity  field  typically  having  velo¬ 
city  excesses  at  inner  radii  and  velocity  defects  at  the  outer  radii. 
This  type  of  velocity  field  has  a  complex  harmonic  content  distri¬ 
bution  and  its  effect  on  blade-rate  noise  cannot  be  predicted 
without  detailed  examination  of  the  wake  and  the  radial  distribu¬ 
tion  of  propeller  blade  forces. 

Techniques  are  available  for  computing  unsteady 
forces  and  skew  distributions,  and  these  methods  range  from 
low-aspect  ratio  approximations  to  unsteady  airfoil  theory  to 
complete  unsteady,  lifting-surface  methods  However,  since  no 
method  was  available  to  systematically  determine  an  optimum 
skew  distribution  for  the  reduction  of  unsteady  forces,  the  pro¬ 
peller  skew  optimization  program  SKEWOPT  (Greenblatt,  1978; 
and  Parsons  and  Greenblatt,  1978)  was  developed.  SKEWOPT 
determines  a  quadratic  or  cubic  skew  distribution  using  an 
optimization  technique  which  finds  the  set  of  parameters  for 
which  a  user-defined  linear  combination  of  the  unsteady  force 
and  moment  amplitudes  are  minimized.  SKEWOPT  was  written 
for  use  in  ship  propeller  design;  however,  the  intention  was  that, 
alter  minor  modifications  such  as  the  inclusion  of  higher  order 
harmonic  groups,  it  would  be  suitable  for  torpedo  and  submarine 
propeller  design.  Since  the  force  calculation  method  in 
SKEWOPT  was  not  sufficiently  documented,  the  method  was 
replaced,  and,  due  to  the  difficulties  encountered  in  modifying 
SKEWOPT,  a  new  program  (SKEWOPT-2),  patterned  after 
SKEWOPT,  was  written. 

In  the  following  discussion,  the  method  used  in 
SKEWOPT-2  for  determining  the  optimum  skew  distribution  for 
a  propeller  operating  on  an  axisymmetric  body  with  upsteam  con¬ 
trol  surfaces  will  be  described.  The  method  optimizes  the  skew 
(combination  of  warp  and  rake)  such  that  the  unsteady  propeller 
blade-rate  forces  are  minimized.  The  forward  propeller  of  a 
counterrotating  propeller  set  win  be  used  to  demonstrate  use  of 
the  method.  While  the  discussion  and  results  presented  deal 
specifically  with  a  marine  propeller  application,  it  is  felt  that  the 
method  is  applicable  to  any  propeller  operating  in  a  non-uniform 
inflow  velocity  field. 


The  velocity  data  used  in  this  study  was  obtained 
from  wind  tunnel  tests  conducted  in  the  Northrop  (NORaIR) 
subsonic  wind  tunnel  located  in  Hawthorne,  California  (see  Nel¬ 
son  and  Fogarty,  1977).  In  order  to  obtain  body  drag  measure¬ 
ments,  the  0.6  scale  model  was  strut  mounted  in  such  a  manner  to 
reduce  horizontal  buoyancy  effects,  and  the  tests  were  conducted 
over  a  Reynolds  number  range  of  1.3-2.4x10*  per  foot  The 
boundary  layer  measurements  were  made  using  a  pair  of  pitot 
tube  rakes.  One  rake  contained  four  static  pressure  tubes  and  the 
second  rake  had  eight  total  head  tubes.  To  avoid  the  effects  of 
die  strut,  wake  measurements  were  made  over  an  arc  slightly 
larger  than  90°  on  the  upper  side  of  the  body,  and  the  center  of 
the  measurement  arc  coincides  with  the  fin's  trailing  edge  (body 
had  four  identical  fins). 

The  measured  wake  data,  shown  in  figure  1,  exhibits 
the  velocity  excesses  at  the  inner  radii  and  the  velocity  deficits  at 
the  outer  radii,  produced  by  both  viscous  and  potential  effects, 
typically  found  behind  an  appendage  on  an  axisymmetric  body. 
The  velocity  excesses  are  due  to  the  horseshoe  vortex  formed  at 
the  appendage /body  intersection,  and  the  velocity  deficits  at  the 
outer  radii  are  due  to  the  fin's  tip  vortex  (see,  for  example,  Gree¬ 
ley  and  Milewski,  1986). 

Velocity  fields  of  this  type  play  an  important  part  in 
the  design  of  wake-adapted  propellers  because  the  spatial  varia¬ 
tions  result  in  the  generation  of  unsteady  forces.  The  velocity 
data  in  figure  I  represents  a  time  averaged,  spatially  varying 
field,  and  since  the  spatial  velocity  distribution  is  periodic  and 
continuous,  it  may  be  represented  in  terms  of  a  Fourier  series. 
For  example,  the  axial  component  of  the  velocity,  at  a  position 
(M),  can  be  expressed  as 


-  — p-  ♦  £  [a„(r)cos(m0)  ♦  bjr)  sin(rro)j 
ajr)-ibm(r)  exprimi?)  j 

J  J 

where  Ref  )  denotes  the  real  part,  and  the  Fourier  coefficients 
are  defined  by 

•*(r>-7/V*d» 

*  -4P  « 

ajr)  -  ±/  ^-cosfmrf)  d#  (2) 

Mr)  -  ^/-^-siiKnrf)  <U 

Since,  the  term  a,(r)  does  not  vary  in  t,  it  is  associated  with  the 
steady  state  thrust  and  torque,  and  the  additional  terms  are 
sinusoidal  fluctuations  of  the  inflow  velocity  which  produce  the 
unsteady  forces  and  moments. 
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THE  SKIW  OPTIMIZATION  PROBLEM 
VELOCITY  FIELD 

In  the  design  of  wake-adapted  propellers,  it  is  impor¬ 
tant  that  the  inflow  velocity  distribution  be  properly  specified. 
The  state-of-the-art  in  boundary  layer  computations  is  such  that 
the  inflow  velocity  field,  for  a  body  having  appendages  located 
npwream  of  the  propulsor,  should  be  determined  experimentally. 
Even  though  circumferentially  averaged  velocity  profiles  are  suf¬ 
ficient  for  propeller  design  calculations,  the  calculation  of 
unsteady  forces  requires  that  both  the  radial  and  circumferential 
distributions  of  the  wake  be  considered. 


PROPELLER  GEOMETRY 

The  propeller  geometry  chosen  for  the  present  study  is 
that  for  the  forward  propeller  of  a  counterrotating  propeller  set 
designed  using  the  method  developed  by  Nelson  (1972.  1975). 
The  design  utilized  the  circumferential  mean,  inflow  velocity 
field  measured  by  Nelson  and  Fogarty  (1977)  and  the  parameters 
given  in  table  I.  The  propeller  geometry,  shown  schematically  in 
figure  2,  was  determined  using  the  lifting-line  portion  of  Nelson's 
design  method.  The  calculated  performance  parameters  for  the 
counterrotating  propeller  set  are  given  in  table  I  while  the  details 
of  the  forward  propeller  geometry  ere  given  in  table  2. 
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Once  the  velocity  Add  and  propeller  geometry  are 
determined,  the  problem  of  determining  an  ’optimum*  skew  dis¬ 
tribution  requires  the  formulation  of  a  nonlinear  programming 
problem  which  includes  an  unsteady  force  calculation  method. 
First,  the  calculation  of  the  unsteady  forces  will  be  considered. 


UNSTEADY  FORCE  CALCULATION  METHOD 
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During  th<?  skew  optimization  process,  unsteady  forces 
will  be  calculated  many  times,  therefore,  it  would  be  desirable  to 
use  an  efficient,  inexpensive  method.  The  original  version  of 
SKEWOPT  (Grccnblatt,  1978)  had  both  a  two-dimensional, 
unsteady  and  a  more  time  consuming  lifting-line  method  available 
to  calculate  the  forces.  To  overcome  both  documentation  and 
computation  problems,  the  original  SKEWOPT  force  calculation 
methods  were  replaced  by  a  method  developed  by  Thompson 
(1976).  His  method  divides  the  propeller  blade  into  strips  which 
are  considered  two-dimensional  airfoils.  Included  in  the  method 
are:  a)  the  two-dimensional  unsteady  airfoil  theories  of  Scars 
(1941)  and  Horlock  (1968)  which  allow  consideration  of 
sinusoidal  velocity  fluctuations  normal  and  parallel  to  the  inflow 
velocity;  and  b)  corrections  to  the  blade  lift  force  due  to  the 
presence  of  adjacent  propeller  blades.  The  method  has  been 
extended  to  include  the  effects  of  camber  using  the  method  of 
Naumann  and  Ych  (1973)  and  to  calculate  the  total  force  and 
moment  on  the  propeller.  A  brief  description  of  the  method  fol¬ 
lows  (for  details  of  the  method  see:  Thompson,  1976;  Mautner, 
1987b;  and  Mautner  and  Blaisdeil.  1987) 

Expressions  for  the  unsteady  lift  acting  on  a  two- 
dimensional  airfoil  due  to  a  periodically  varying  free  stream  velo¬ 
city  were  developed  by  Scars  (1941)  and  Horlock  (1968).  They 
assumed  that  the  airfoil  can  be  represented  by  an  isolated  flat 
plate,  and  that  the  flow  is  incompressible,  inviscid  and  irrota- 
tional  except  for  surfaces  of  distributed  singularities.  The  airfoil 
is  replaced  by  a  distribution  of  vortices,  and  the  shed  vorticity  is 
assumed  to  lie  on  a  line  parallel  to  the  unsteady  free  stream  velo¬ 
city.  Also,  the  inflow  velocity  and  the  velocity  induced  by  the 
vortex  system  must  satisfy  the  boundary  condition  that  there  is 
no  flow  through  the  plate.  Their  formulation  results  in  an 
integral  equation  describing  the  vortex  distribution,  and  the  solu¬ 
tion  of  the  integral  equation  yields  the  total  velocity  on  the  plate's 
surface.  Using  Bernoulli's  equation,  the  pressure  distribution  on 
the  upper  and  lower  surfaces  of  the  plate  arc  obtained,  and  the 
unsteady  lift  is  determined  from  the  pressure  difference. 

In  order  to  apply  two-dimensional,  unsteady  airfoil 
theory  to  propellers,  the  propeller  blade  is  divided  into  strips  of 
width  Art  as  indicated  in  figure  2.  Each  (Hade  element  is  then 
treated  as  an  isolated  two-dimensional  airfoil  having  at  its  mid¬ 
point  the  relative  velocity  V,,.  The  velocities  which  determine 
V„  are  the  propeller  rotation  flrj  ,  the  mean  axial  inflow  velocity 
V,(r^)  and  the  down  wash  velocity  (which  for  wake-adapted 
propellers  includes  induced  and  interference  velocities).  The 
resulting  angle  between  Or,  and  V„  defines  the  local  pitch  angle 
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The  varying  axial  velocity  field  incident  upon  the 
rotating  blade  row  is  resolved  into  components  normal  and  paral¬ 
lel  to  V„  .  The  assumption  of  a  lightly  loaded  propeller  is  made, 
the  velocity  fluctuations  are  expressed  in  blade  fixed  coordinates 
for  each  n-th  harmonic  component,  and  the  phase  shift  6) 
depends  on  both  the  skew  angle  p,  and  the  blade  spacing 
m-2»l/Nb- 


Adding  the  contributions  of  the  gnst  velocities  to  the 
unsteady  lift-  and  replacing  the  factor  of  2tr  with  the  Made  sec¬ 
tion  lift  curve  slope  to  account  for  the  additional  lift  due  to  adja¬ 
cent  blades,  one  obtains  for  the  lift  and  moment  on  a  Made  sec¬ 
tion 


where  the  reduced  frequency  is  k^  ■  nfi  ^C/V^  .  The  lift  curve 
slope  can  be  determined  experimentally  or  approximated,  for 
example,  by  the  method  of  Weinig  described  in  Wisliccnus  (1965). 

As  the  propeller  rotates,  blade  position  dependent  x 
and  y  components  of  the  lift  vector  and  the  moment  arm  are  gen¬ 
erated.  Therefore,  in  addition  to  the  thrust  and  torque,  there  are 
x  and  y  force  and  moment  components.  Also,  it  has  been  shown 
that  the  only  harmonics  contributing  to  the  thrust  and  torque  are 
those  of  order  mNb  (some  multiple  of  blade  number),  and  the 
only  harmonics  contributing  to  the  side  forces  and  moments  are 
of  order  mN^AI.  Since  the  contributions  from  the  different  har¬ 
monic  groups  oscillate  at  different  frequencies  they  should  be 
considered  separately  and  doing  so  results  in  the  following 
expressions  for  the  unsteady  forces  and  moments  due  to  a  partic¬ 
ular  harmonic  group  m. 

(5) 

jil  I  4  1 

T.w-Nb  £  f Mu-uj  I ,  <6* 

Fj**'  -  *NbEj  je  (8) 

T«  -  i  £  {  (Mil-vo),'  (»W.l)t  <’> 

Tj"*  ■  laNbsj  j*  <10' 

In  addition  to  the  above  forces  and  moments. 
Thompson's  method  was  extended  (Mautner.  1917b;  Mautner 
and  Blaisdeil,  1987)  to  include  calculation  of  the  maximum  side 
force  and  bending  moment.  Briefly,  the  forces  and  moments  are 
expressed  in  complex  form,  and,  via  a  coordinate  rotation,  the 
derived  equations  for  the  total  side  force  and  bending  moment 
vector  trace  an  ellipse  where  the  semi-major  axis  yields  the  max¬ 
imum  amplitude  of  the  side  force  and  bending  moment  The 
equations  are 


on 


(12) 
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The  subscripts  1  snd  2  used  in  the  force  and  moment  equations 
(5)-,  t  indicate  the  following  forms  of  the  lift  and  moment  equa¬ 
tions 


optimization  search  to  be  carried  out  in  either  a  two-  or  three- 
dimensional  space. 
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where  the  index  n  takes  on  the  values  mN*  .  mNbtl  or  mNb-l. 
This  completes  the  formulation  of  unsteady  force  and  moment 
equations.  A  description  of  the  skew  optimization  method  is 
given  in  the  next  section. 


SKEW  OPTIMIZATION  MODEL 

To  determine  the  optimum  skew  distribution  the 
above  force  calculation  method  was  incorporated  into  a  nonlinear 
programming  problem.  Due  to  the  fact  that,  in  general,  all  forces 
cannot  be  minimized  simultaneously,  a  scalar  cost  function 
formed  from  the  weighted,  linear  combination  of  the  forces  and 
moments  is  minimized.  The  cost  function  Fc  is 


Wj"*F.u 


.05  F. 


.05  F, 


where  the  weights.  W/”1.  are  normalized  such  that  their  sum  over 

both  m  and  j  (-1, _ ,4)  equals  1.  The  normalization  factors  (i.c. 

the  steady  state  thrust  and  torque.  and  1",  )  arc  arbitrary;  how¬ 
ever.  the  WjM  are  chosen  to  place  emphasis  on  the  suppression  of 
particular  forces  or  moments. 

Since  Fc  depends  upon  the  skew  distribution,  d(r),  one 
could  solve  for  the  optimum  skew  distribution  via  a  variational  or 
n-dimensional  parameter  optimization  technique.  However,  a 
more  feasible  approach,  and  the  one  used  here,  is  to  use  a  few 
parameters  in  describing  the  skew  distribution  and  perform  an 
optimization  search  in  a  limited  parameter  space.  To  accomplish 
this,  the  skew  distribution  is  represented  by  either  s  cubic  or 
quadratic  distribution  having  a  straight  line  section  with  *"0. 
The  skew  distribution  (see  Gcenblatt,  1971;  Mautner  and  Blais- 
dell,  1987),  illustrated  in  figure  3,  is  described  by 


The  propeller  skew  distribution  problem  has  now  been 
formulated  as  an  optimization  problem  in  terms  of  a  few 
geometric  parameters.  In  order  to  obtain  a  feasible  propeller 
geometry,  it  is  necessary  to  place  some  restrictions,  such  as  a 
maximum  allowable  tip  skew,  on  the  geometric  parameters.  In 
doing  so,  the  design  problem  becomes  a  constrained,  nonlinear 
optimization  problem  where  the  optimization  search  is  restricted 
to  finding  the  set  of  parameters  which  minimizes  the  cost  func¬ 
tion  Fe  while  satisfying  all  of  the  constraints  placed  on  the  prob¬ 
lem.  The  constraints  used  are  given  in  table  4  and  are  checked 
for  violation  at  a  given  point  in  the  parameter  space  using  the 
method  developed  by  Creenblatt  (1978).  The  additional  con¬ 
straint  that  the  propeller  blade  should  not  curve  forward  has 
been  added. 

The  constrained,  nonlinear  optimization  problem  is 
represented  by 


min  F,(X)  subject  to  Gj  (X)  >  0  (17) 


where  Xiti  vector  in  parameter  space  which  determines  the 
skew  distribution  and  Gj  (X.)  describes  the  constraints.  There  ar: 
many  techniques  for  solving  the  unconstrained  minimization 
problem  (see  Parsons,  1973);  however,  only  a  few  methods  attack 
the  constrained  problem  directly.  One  useful  technique  is  to  con¬ 
vert  the  constrained  problem  into  an  unconstrained  problem  and 
then  use  an  unconstrained  optimization  method.  This  can  be 
accomplished  with  the  use  of  an  external  penalty  function  which 
is  added  to  a  cost  function  whenever  a  constraint  is  violated  (i.c. 
Gj(5£)  <  0  ).  The  penalty  function  is 

Pe(X,Ak)  -  F,<3£)  -  Ak  £  min  [Gj(J£).  01  (IS) 

i 

An  unconstrained  optimization  technique  can  be  applied  to 
Pe(X  ,Ak).  If  no  constraint  is  violated  no  penalty  is  added  and 
the  penalty  function  is  the  same  as  die  cost  function.  Since  the 
penalty  added  is  proportional  to  the  constraint  violation,  the 
optimisation  method  should  be  forced  towards  a  feasible  region 
where  no  constraints  are  violated.  This  will  be  the  case  as  long 
as  the  multipiieative  factor  Ak  is  large  enough  (*1024);  however, 
if  Ak  it  too  small,  the  search  may  tend  toward  an  infeasible 
region.  The  optimization  technique  chosen  for  use  in  this  method 
is  the  Ncldcr-Mcad  simplex  search  method  (Ncldcr  and  Mead. 
1963),  and  a  detailed  description  of  the  method  as  coded  in  the 
current  computer  program  is  given  by  Mautner  and  Blaisdcll 
(IM7). 


tHr)»Alr*+Aui**Aar+A4  or  A^A^vA,  r.JrSr, 
p(r)-0  rkSrSr,(16) 

The  skew  distribution  given  above  has  five  free 
parameters  (  A, ,  A, ,  A, .  A« ,  r, );  however,  it  is  more  meaning- 
Ibl  to  the  propeller  designer  to  use  parameters  that  have  physical 
meaning  instead  of  using  these  polynomial  coefficients.  The 
parameters  chosen  for  use  in  the  current  method  are  the  skew  at 
the  propeller  tip  the  muting  skew  slope  S,  -  v'(r.),  the  start- 
tag  radius  r„  and  the  skew  dope  at  the  tip  &.  Additionally,  the 
physical  restriction  that  9(0*0  j,  made  so  that  the  skew  distri¬ 
bution  is  continuous  thereby  reducing  the  number  of  free  param¬ 
eters  by  one.  Another  necessary  restriction  is  that  the  skew  dis¬ 
tribution  be  smooth  which  implies  that  S,  ■  B9(r,)/dr-0  if  r,  >  rk. 
r  there  is  no  straight  Itaa  taction  r,  -  S,  is  not  restricted  but  r, 
is  fixed.  In  either  earn,  the  eumber  of  free  parameters  is  reduced 
» three  for  a  cubic  distribution  and  two  Ibr  a  quadratic  distribu¬ 
tion  and  results  in  the  four  possible  skew  distribution  models 
presented  ta  table  1.  Specification  of  diass  models  alewt  the 


The  computer  program  which  solves  the  above  non¬ 
linear  programming  problem  is  an  interactive  program  intended 
for  routine  propeller  design  work.  The  program  has  an  interactive 
input  method  to  accept  propeller  data  and  optimization  parame¬ 
ters,  and  the  program  can  be  restarted,  for  example,  with  dif¬ 
ferent  constraint  values.  In  addition  to  the  optimization  mode,  a 
test  mode  is  available  for  the  calculation  of  unsteady  forces  for  a 
given  skew  distribution. 


NUMERICAL  EXAMPLES 
FOUblER  ANALYSIS 

In  addition  »  the  specification  of  propeller  geometry, 
a  preliminary  step  in  the  calculation  of  unsteady  forces  and  skew 
distributions  is  an  analysis  of  dm  input  wake.  Examination  of  die 
annaailj  fores  aquations  shows  that  a  Fourier  analysis  of  the 
wake  la  required,  snd  that  only  certain  harmonics  of  the  input 
wake  will  coutrfoute  to  the  unsteady  forces  and  momenta.  Tha 
Fourtar  analysis  ef  the  currant  four  cycle  wake  shows  dm  dta 


harmonic  numbers  of  interest  are  4,  S,  12,  16, _ (all  others  have 

zero  magnitude).  Also,  since  a  six  bladcd  propeller  is  being  used 
in  this  example,  the  harmonics  of  concern  include  6,  12.  18,  24. 

_  Thus  for  a  six  bladcd  propeller  operating  in  a  four  cycle 

wake,  the  only  harmonics  of  interest  are  12,  24,  36, _ (If,  for 

example,  a  seven  bladed  propeller  was  specified,  one  would  con¬ 
sider  side  forces  due  to  harmonics  of  order  mN^±  I.)  The  radial 
distribution  of  the  Fourier  coefficients  and  phase  angles  for  the 
12th  and  24th  harmonics  are  given  in  figures  4,  S  and  6.  The 
results  show  both  the  dominance  of  the  12th  harmonic,  especially 
in  the  region  of  r/R0  <  0.3,  and  the  rapid  approach  to  a  nearly 
zero  magnitude  of  the  24th  and  higher  harmonics.  For  example, 
the  peak  value  of  c*  for  the  24th  harmonic  is  2.3  times  lower  than 
the  peak  c.  for  the  1 2th  harmonic.  The  phase  angle  distributions 
(figure  6)  show  that  the  phase  for  the  24th  harmonic  is  nearly 
constant  over  radius  while  the  phase  angle  of  the  12th  harmonic 
has  a  nearly  constant  positive  magnitude  for  0  <  r/Ra  <  a  0.3,  it 
undergoes  a  sharp  phase  shift  at  r/RB »  0.3,  and  then  has  a 
nearly  constant  negative  magnitude  over  a  0.6  <  r/R3  <  a  0.9. 
Finally,  since  the  harmonic  numbers  of  concern  are  12,  24,  36, 
—  the  only  forces  that  need  to  be  calculated  are  the  unsteady 
thrust  and  torque. 


UNSTEADY  FORCES 

To  illustrate  use  of  the  unsteady  force  calculation 
method,  the  six  bladed  propeller  geometry  and  operating  parame¬ 
ters  (tables  1  and  2)  along  with  the  12th  and  24th  harmonic  dis¬ 
tributions  were  provided  as  inputs  to  the  method.  The  calculated 
ndia!  distributions  of  F.  and  T.  are  given  in  figures  7  and  8  for 
the  12th  and  24th  harmonics.  The  radial  distributions  for 
mNb»  12  show  that  regions  of  large  forces  occur  in  both  the  inner 
and  outer  portions  of  the  propeller  blade.  Also,  there  is  a  distinct 
minimum  force  region,  located  at  r/R  *  0.68,  which  coincides 
with  the  minimum  velocity  defect/excess  region  of  r/  Rg  *  0.33 
shown  in  figure  1.  The  forces  for  the  24th  harmonic  are  signifi¬ 
cantly  lower  that  those  for  1 2th  harmonic  and  increase  in  magni¬ 
tude  with  increasing  radius.  This  result  follows  directly  from  the 
lower  harmonic  magnitudes  for  mNb«24.  Radial  distributions  of 
forces  such  as  that  shown  in  figures  7  and  8  provide  information 
relevant  to  what  portion  of  the  wake  is  most  important  in  the 
production  of  unsteady  forces  and  in  the  possible  reduction  of 
these  forces  by  appropriate  appendage/vehicle  design.  A  detailed 
comparison  of  the  wakes  behind  various  appendage  geometries 
and  the  resulting  unsteady  forces  can  be  found  in  Maumer 
(1987a). 


SKEW  DISTRIBUTIONS 

To  demonstrate  the  types  of  skew  distributions  which 
can  be  obtained  using  the  current  optimization  method,  the  four 
models  summarized  in  tables  3  and  4  were  used  to  calculate 
'optimum'  skew  distributions  and  the  magnitude  of  the  unsteady 
thrust  and  torque  associated  with  the  skew  distributions.  To  pro¬ 
vide  a  reference,  the  magnitude  of  the  unsteady  thrust  and 
torque  for  the  unskewed  propeller  was  calculated.  Next,  each  of 
the  four  skew  distribution  models,  with  essentially  no  limits  on 
the  constraints,  was  specified.  For  example,  the  limit  on  max¬ 
imum  tip  skew  was  2*/Nb.  The  calculated  unsteady  thrust  and 
torques  are  given  in  able  3  where  the  reduction  in  forces 
obtained  for  each  skew  distribution  is  represented  by  FV  F,  and 
T„  T,  where  the  bar  represents  the  wtskewtd  values.  The  calcu¬ 
lated  skew  distributions  arc  plotted  in  figures  9-11  where  the 
letters  correspond  to  the  particular  cases  listed  in  table  3. 

It  is  dear  from  an  examination  of  the  calculated  forces 
given  in  table  3  (cases  A,  i,  C.  F)  that  a  significant  reduction  of 
the  total  unsteady  thrust  and  torque  was  achieved  for  the  24th 
harmonic  irregardkss  of  the  skew  distribudoe  model  wed.  How¬ 
ever,  Car  the  12a  harmonic,  the  quadratic  distributions  *l  and 


*2  (A,B)  and  the  cubic  *4  (F)  distribution  produced  signifi¬ 
cantly  less  reduction  in  the  magnitude  of  both  the  thrust  and 
torque.  The  largest  force  reduction,  for  mN,,-l2.  occurs  for  cubic 
*3  (C)  which  also  has  a  shape  (figure  9)  that  deviates  drastically 
from  the  general  shape  of  the  other  three  skew  distributions 
(A,B,F).  The  resulting  shape  of  cubic  «*3  (C)  points  out  a  poten¬ 
tial  problem  with  the  specification  of  unlimited  constraints.  While 
the  method  attempts  to  provide  the  skew  for  minimum  force 
reduction,  it  may  provide  a  skew  distribution  which  is  not  feasi¬ 
ble  in  terms  of  final  propeller  geometry.  The  results  of  applying 
constraints  to  the  cubic  »3  model  arc  shown  in  figure  10.  When 
the  starting  skew  slope  is  limited,  the  resulting  skew  distribution 
and  force  reduction  fall  in  tine  with  the  results  obtained  for  the 
other  models  (figure  9). 

As  shown  by  equation  (IS),  the  skew  optimization 
problem  involves  the  linear  combination  of  weighted  forces.  The 
weights  Wj  are  specified  by  the  propeller  designer  and  thus  can 
be  chosen  to  place  more  emphasis  on  the  reduction  of  a  particu¬ 
lar  force  in  a  particular  harmonic  group.* To  illustrate  the  effect  of 
varying  the  weights  W,  the  cubic  #3  and  cubic  »4  models  were 
used.  In  addition  to  the  Wj»l  used  in  the  previous  calculations, 
the  weights  were  increased  to  S  and  10  for  the  12th  harmonic 
while  maintaining  Wt*l  for  the  24th  harmonic.  The  results 
presented  in  able  3  and  figure  1 1  show  very  little  change  in  the 
force  magnitudes  and  the  shape  of  the  skew  distributions.  This 
result  is  not  surprising  due  to  the  dominance  of  the  12th  har¬ 
monic  magnitude  (figures  4  and  5). 

As  mentioned  previously,  the  original  motivation  for 
the  introduction  of  skew  into  marine  propeller  design  was  the 
reduction  of  both  blade-rate  noise  and  vehicle  vibration.  The 
basic  idea  mi  to  introduce  a  phase  shift  in  the  local  periodic  lift 
forces  over  the  radial  extent  of  the  blade.  The  phase  shift  should 
be  specified  such  that  the  local  forces,  at  the  inner  and  outer 
radii,  are  acting  in  opposition.  To  illustrate  the  results  of  applying 
the  force  equations  to  obtain  a  phase  shift,  the  radial  distribution 
of  the  real  an  imaginary  parts  of  unsteady  thrust  distributions  are 
plotted  in  figures  12  and  13  for  mN,,- 12.  The  figures  compare 
the  radial  distribution  of  unsteady  thrust  for  the  case  of  no  skew 
with  the  force  distribution  obtained  using  the  cubic  *3  (C)  skew 
distribution  model  The  radial  distribution  of  both  the  real  and 
imaginary  para  of  F,  demonstrate  the  acquisition  of  a  phase  shift 
in  the  unsteady  lift  by  the  targe  changes  in  magnitude  of  the 
unsteady  force.  Due  to  die  character  of  the  12th  harmonic  distri¬ 
bution,  it  is  not  possible  to  obtain  a  phase  shift  in  the  real  part  of 
F,  such  that  the  forces  at  the  inner  and  outer  radii  are  acting  in 
opposition;  however,  the  forces  are  acting  in  opposition  for  the 
the  imaginary  part  of  Fr  The  results  in  figures  12  and  13  show 
that  the  forces  are  shifted  from  their  original  condition  thus  pro¬ 
ducing.  at  least,  some  unsteady  force  reduction.  Similar  results 
are  obtained  for  the  unsteady  thrust  associated  with  mN^-24  and 
the  ensteady  torque  (not  shown  here)  for  both  the  12a  and  24a 
harmonics. 

The  close  similarity  in  most  of  the  skew  distributions 
(figure  9-11)  indicates  aat  satisfactory,  but  not  necessarily 
optimum,  propeller  geometries  can  be  obtained  without  having  to 
place  severe  restrictions  on  the  skew  distribution  model  con 
strain  ts.  However,  the  magnitudes  of  the  unsteady  forces  and 
moments  listed  in  table  3  demonstrate  aat  for  an  actual  propeller 
design  the  effect  of  parameters  such  as  skew  slope  and  tip  skew 
need  to  be  examined  closely  for  its  effect  on  each  harmonic 
group. 
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SUMMARY 

A  propeller  skew  optimization  program,  based  upon 
Green  Matt's  method,  has  been  developed.  The  method  uses  a 
quadratic  or  cubic  skew  distribution  model,  and  a  parameter 
search,  subject  to  constraints,  is  performed  to  determine  the 
skew  distribution  which  minimizes  a  linear  combination  of  the 
unsteady  forces  and  moments.  It  should  be  remembered  that  this 
program  requires  user  judgement  in  the  specification  of  the  skew 
distribution  model  and  constraints. 

The  two-dimensional,  unsteady  airfoil  theory  used  in 
the  unsteady  force  calculation  method  may  not  be  as  accurate  as 
three-dimensional,  unsteady  lifting-surface  methods  in  predicting 
unsteady  loads  (Le.  for  low  aspect  ratio  blades).  However,  if  the 
errors  are  consistent  so  that  the  method  predicts  the  correct 
trends  in  the  unsteady  forces,  the  resulting  skew  distribution 
should  be  reasonably  accurate.  It  should  also  be  noted  that  the 
original  intention  of  incorporating  Thompson’s  (1976)  force  cal¬ 
culation  method  into  the  optimization  program  was  to  provide  a 
well  documented  and  efficient  method  for  use,  at  least,  during 
the  development  phase  of  the  skew  optimization  program. 
Although.  Thompson  (1976)  obtained  good  agreement  between 
theory  and  experiment  and  the  current  method  provides  reason¬ 
able  skew  distributions,  the  need  exists  for  additional  program 
enhancements. 

One  problem  inherent  in  the  application  of  skew 
involves  the  use  of  wake  data  obtained  without  a  propeller 
present  While  the  propeller  design  method  of  Nelson  (1972, 
197$)  calculates  changes  in  the  circumferential  mean,  inflow 
velocity  field  due  to  presence  of  the  propulsor  (Le.  induced  and 
interference  velocities),  the  uncorrected,  spatially  varying  inflow 
is  used  to  determine  the  unsteady  forces.  It  is  known  that  the 
presence  of  a  propulsor  will  cause  a  change  in  the  streamlines 
due  to  acceleration  of  the  flow,  that  there  may  be  additional 
unsteadiness  due  to  the  relative  motion  of  the  blades  and  that  the 
propeller  will  change  the  amplitude  and  phase  of  the  incident 
distortion  velocity.  From  these  few  facts  it  is  apparent  that  the 
incident  wake  should  be  corrected  for  propulsor  effects. 

Attempts  are  being  made  to  develop  unsteady  force 
calculation  methods  which  account  for  the  effects  of  both  the 
presence  of  the  propulsor  and  the  Made  skew  on  the  harmonic 
content  of  the  incident  wake.  One  such  effort  is  that  by  Zierke 
(19SS).  He  has  pointed  out  that  both  the  amplitude  and  phase 
angle  of  the  inlet  distortion  will  be  modified  by  the  blade  skew 
(warp  affects  both  amplitude  and  phase  and  rake  affects  the 
amplitude),  and  that  skew  is  a  measure  of  when  a  blade  section 
first  encounters  a  disturbance.  Zierke  has  developed  an  unsteady 
force  calculation  method  which  both  accounts  for  the  changes  in 
amplitude  and  phase  of  a  disturbance  due  to  the  presence  of  a 
rotating  blade  and  uses  leading  edge  skew  rather  than  mid-chord 
skew  to  describe  the  Made.  An  obvious  improvement  to  the 
current  method  would  be  the  incorporation  these  kinds  of 
corrections. 
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Table  1.  Theoretical  Performance  of  the 
Counterrotating  Propeller  Set 


Propeller  radius  (in) 

R,-R 

8.: 

RPM 

RPM 

1400 

Propulsive  coefficient 

PC 

0.92S 

Thrust  deduction  factor 

T 

0.141 

Vehicle  speed  (knots) 

V. 

43 

Drag  coefficient 

Cp 

0.113 

Advance  ratio 

J 

2.12 

Blade  surface  cavitation 

a 

0.75 

Thrust  coefficient 

Ct 

0.21S 

Torque  coefficient 

Cq 

0.137 

Blade  tip  speed 

v*/v.v 

1.731 

V«p/V.)A 

1.604 

Number  of  blades 

Nb.r  *  Nb.A 

6x4 

Hub  radius  (in) 

rk 

3.44 

Body  radius  (in) 

R> 

10.5 

Subscripts  *  F  *  Forward  and  A  -  After  Propeller 


Table  2.  Propeller  Geometry  and  Operating  Characteristics 


Table  3.  Skew  Distribution  Models. 


Model 

m 

Free 

Parameters 

Fixed 

Parameters 

1 

Quadratic 

M 

mm 

2 

Quadratic 

3 

Cubic 

4 

Cubic 

*«  .  s, ,  r. 

UJI 

Table  4.  Skew  Distribution  Constraints 


HOKfMAf  NSIONAI  VI  l  OOlt  Vli  »l 


Table  S.  Calculated  Unsteady  Force  Reduction 
and  Tip  Skew. 


Skew  Model 
(see  Table  3) 

mN, 

F.* 

-12 

t; 

mNt 

F. 

-24 

T. 

l 

u  » 
% 

wR 

No  Skew 

38.1 

12.8 

13.4 

5.4 

0.0 

1 

F./?.* 

VT* 

F„F, 

Ti/T, 

No  Skew 

1.00 

1.00 

1. 00 

1.00 

0.0 

1 

(A)  Quadratic  *1 

0.70 

0.55 

0.17 

0.07 

59.6 

I 

(B)  Quadratic  #2 

0.65 

0.52 

0.25 

0.17 

57.4 

l 

(C)  Cubic  *3 

0.09 

0.19 

0.09 

0.07 

58.6 

1 

(D)  Cubic  *3+ 

0.31 

0.21 

0.14 

0.11 

59.7 

1 

(E)  Cubic  *3++ 

0.42 

0.29 

0.15 

0.11 

59.1 

1 

(F)  Cubic  *4 

0.50 

0.37 

0.13 

0.07 

59.8 

1 

(G)  Cubic  «3 

0.08 

0.02 

0.11 

0.09 

59.5 

5 

(H)  Cubic  *3 

0.08 

0.02 

0.11 

0.09 

59.5 

10 

(I)  Cubic  *»4 

0.39 

0.26 

0.12 

0.09 

59.6 

5 

(J)  Cubic  -4 

0.42 

0.29 

0.14 

0.09 

59.2 

10 

Constraints 

Tip  Skew  -  Max  -  60*  and  Min  -  -30° 

Maximum/  Minimum  Skew  Slope  -  ±5.38 
Start  Radius  •  rfc-3.44  in 

Notes 

•  Units  (F)orce  (thrust)-! bs;  (T)orque-fl-lbs 
*•  Mid-chord  tip  skew  -  deg 

#  Bar  indicates  no  skew  thrust  and  torque 
+  Starting  slope  constraint  -  -I  <  S^  £  +1 

++  Starting  slope  constraint  •  -0.25  S  S,  5  +025 
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FOURIER  COEFFICIENTS  b„.e„  "d  nj  FOURIER  COEFFICIENTS  a 


'  °r  I  r  •=  r( 


#  -0. 8  ■  \ 

<  mN0  =  24  \ 

L  - *  V."JIH 

-2.4  - 


T—  \  S  S  ^1  -4.0' - 1 - 1 - ' - 1 - 1 - 1 - 1 - 1 - 1 - 1 

/  .  *  df  1  r  =  rt  0  20  0.36  0.52  0.68  0  84  V00 

f-N.  RADIAL  POSITION.  r/R. 


Figure  3.  Definition  of  the  Skew  Distribution  Model  Parameters. 


Figure  6.  Variation  of  the  Fourier  Coefficient  Phase  Angle  with 
Radial  Position  for  the  12th  and  24th  Harmonics. 


mNs  =  12  F-’ 

\ 

/ 

F:~J 

\  / - n 

1 

\  / 

1 

A 

\  / 

T.-y" 

\  /  — 

o  20  0  36  0  62  0  68  0  84  1  X) 

RAOIAL  LOCATION  r  R. 


6  0.52  0.68 


NON-DIMENSIONAL  RADIAL  POSITION.  r/R 

Figure  7.  Variation  of  the  Unsteady  Thrust  and  Torque  with 
Radial  Position  for  the  12th  Harmonic. 


gure  4.  Variation  of  Fourier  Coefficient  Amplitude  With  Radial 
Kition  for  the  12th  Harmonic. 


-0041— 
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0.20  036 


052  068  0  84  1  00 

RAOIAL  LOCATION  r.  R. 


0  36  0.52  0  68  0  84  1  00 

NON-DIMENSIONAL  RAOIAL  POSITION.  r/R 


Figure  S.  Variation  of  Fourier  Coefficient  Amplitude  With  Radial 
Position  for  the  24th  Harmonic. 
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Figure  8.  Variation  of  the  Unsteady  Thrust  and  Torque  with 
Radial  Position  for  the  24th  Harmonic. 


MID-CHORD  SKEW  ANGLE  (deg) 


Figure  9.  Calculated  Mid-Chord  Skew  Distributions. 


MIO-CHORO  SKEW  ANGLE  (deg) 


Figure  10.  Calculacsd  Mid-Chord  Skew  Distribution]  Using  the 
Cubic  * 3  Model  With  Various  Constraint! 


MID-CHORD  SKEW  ANGLE  (deg) 


Figure  II.  Calcuated  Mid-Chord  Skew  Distributions  Using  the 
Cubic  #3  and  *4  Models  for  Various  Values  or  the  Weighting 
Function. 


Figure  12.  Variation  of  the  Real  Part  of  the  Unsteady  Thrust 
with  Radial  Position  for  the  12th  Harmonic. 


RADIAL  POSITION.  r/R 


Figure  13.  Variation  of  the  Imaginary  Part  of  the  Unsteady  Thrust 
with  Radial  Position  for  the  12th  Harmonic. 
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